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Abstract 

We prove weighted-L 00 and pointwise space-time decay estimates for weak solutions of a class of 
wave equations with time-independent potentials and subject to initial data, both of low regularity, 
satisfying given decay bounds at infinity. The rate of their decay depends on the asymptotic 
behaviour of the potential and of the data. The technique is robust enough to treat also more 
regular solutions and provides decay estimates for arbitrary derivatives, provided the potential and 
the data have sufficient regularity, but it is restricted to potentials of bounded strength (such that 
—A — |V| has no negative eigenvalues). 



I. INTRODUCTION 

We study a class of wave equations of the form 

Uu + Vu := d]u - Au + V(x)u = (1) 

where u := u(t,x) : R+ x R 3 — >• R and V(x) is a real potential which does not depend on 
time. We are interested in weak solutions to the initial value problem 

u{0,x) = f(x), d t u(0,x)=g(x) (2) 

with data (/, g) of low regularity satisfying either some weighted-L 00 or pointwise bounds. 
We prove a decay estimate in two versions: in stronger, we show a pointwise decay 

C — 

\u(t,x)\<- ; — -- ; , ,, x - V(t,x)eR+xR 3 (3) 

1 v n ~ (i + |t + |x||)(i + |t- \x\iy- 1 v ; + K J 

with some p > 2 provided the potential V and the initial data /, V/, g are continuous and 
satisfy pointwise bounds 

|y ^(TW k>2 (4) 

with | Vb| < C~l (the value of C p ^ will be specified later) and 

l/WI < ^- — -, \Vf(x)\<- , \g(x)\<-, — m>3. (5) 

\J\ n - ( 1 + | x |)m-i' i (1 + |^|) m ~ (1 + \x\) m KJ 

In the weaker version, all pointwise bounds need to hold only almost everywhere and so we 
prove a weighted-L°° space-time decay estimate 

||(1 + \t+ \x\\)(l + \t- Mlf-M^IUfi^RS) < oo (6) 
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with some p > 2 provided 



V := ||(1 + | ■ |)V|| L oo (R3) < C'l < oo, k > 2 (7) 
(the value of C p ^ will be specified later) and 

fo := 11(1 + I " |) TO -7lU~CR3) < oo, /i := ||(1 + | • \) m Vf\\ L -p>) < oo, 

go := ||(1 + | • |) m ^|U- (R3) < oo, m>3. (8) 

For both kinds of initial data, u G L 1 1 oc (R 1+3 ) represents a weak solution of the wave 
equation ([I]) in the following sense: for any test function (p G C°°(IR 1+3 ) 

dt / d 3 x (□ + V)(p{t, x) u(t, x) = - I d 3 x d t <f(0, x) f(x) + / d 3 x (p(Q, x) g(x). (9) 



The initial data /, V/ and g have just sufficient spatial decay to have finite energy 
and we show by functional analytic methods (sec. IIVI) that this energy stays finite during 
the evolution what guarantees uniqueness of the weak solutions. It follows that the first 
derivatives d t u, Vw G L 2 C L\ oc and hence exist in the weak sense. The second derivatives 
exist only as distributions, although Ou = —Vu G L°° C L\ oc . Therefore, we work in a 
subset of the energy space. 

There is, in principle, no essential problem with rising the regularity of the solution u 
provided /, g and V are more regular and their derivatives have proper decay at infinity. 
Indeed, for (f,g) G C n+1 (M 3 ) x C n (R 3 ) and V G C n (R 3 ) the solution u G C n (M+ +3 ) and for 
n > 2 it becomes classical. In a separate section (sec. IVII) we demonstrate robustness of our 
technique and show how to obtain estimates for derivatives, provided more assumptions on 
the initial data and the potential are given. We also easily reproduce the result of Strauss 
and Tsutaya [l[ in the case n = 2. 

The technique of proving these estimates relies essentially on the Duhamel integral rep- 
resentation formula for solutions of the free wave equation with a source and non-vanishing 
initial data. In order for the Duhamel formula to be applicable we need that the functions 
/, V/ and g are Lebesgue integrable over all or almost all spheres in M 3 . Therefore, we 
first choose continuous data and potential, which guarantee the integrability over spheres 
everywhere. Next, we switch to less regular L°°-spaces where we have integrability over 
spheres only almost everywhere what, with some additional work regarding the measures of 
(null) space and space-time subsets, is still sufficient to complete the proof. It seems to be 
difficult to further generalize the decay theorem by weakening the assumptions and still use 
the same technique of proof based on the Duhamel representation. 

Our motivation for proving these estimates comes mainly from the analysis of nonlinear 
wave equations, where decay proofs rely on the corresponding results for the linearized 
equations. We have been also inspired by several nonrigorous approaches to the problem 
of long-time decay and by numerical observations which consistently show that initial data 
of compact support (or of enough rapid decay) in presence of a potential V(x) ~ \x\~ h 
(\x\ 3> 1) evolve to late-time tails of the form u(t,x) ~ t~ k for big t and fixed x. This has 
been first explained by Ching et al. [2], 0| who approximated the Green's function or later by 
Hod 0] who used some series expansion, while both arguments were non-rigorous. Strauss 
and Tsutaya gave in [l| a rigorous proof of the decay estimate for classical solutions, i.e. 
which are twice continuously differentiable, using the Duhamel representation formula. We 
found the regularity conditions not crucial for controlling the decay of the solution itself and 
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dropped the differentiability conditions for both the data and the potential, thus extending 
the technique to weak solutions. By a careful analysis of well-posedness and uniqueness of 
solutions in the weighted L°° space, involving proofs of measurability and local integrability 
combined with functional analytic approach to uniqueness, we were able to prove the same 
decay estimates (without or with control of the derivatives) and energy conservation for 
rough solutions. This setting covers various cases of cut-off potentials or cut-off initial data 
often discussed in the context of specific approximations or in numerical calculations which 
were not covered by the existing literature. 



Strategy of proof 

We introduce the following notation for solutions of the wave equations. Let Iy be a 
linear map from the space of initial data to the space of solutions of the wave equation 
([I])-®! so that u = Iy(f,g). For wave equations with a source term and null initial data 

n U + Vu = F, u{0,x) = 0, d t u{0,x)=0, (10) 

let's denote the solutions by u = Ly(F), where Ly is a linear map from the space of source 
functions to the space of solutions to the above problem. Due to linearity, the solution u 
of a wave equation with source F and non- vanishing initial data /, g is a sum of these two 
contributions 

U = Ly(F)+Iy(f,g). (11) 

The method we use for proving decay of solutions in presence of a potential is to treat the 
potential term as a source term. If we put the term Vu in flTJ on the r.h.s 

Uu = -Vu, (12) 

we get a pseudo-free wave equation with a source term F = —Vu. Hence, the solution must 
satisfy 

u = I v (f,g) = Lq(-Vu) + I (f,g), (13) 

where the difficulty is that u appears on both sides. By introducing a weighted space-time 
norm II (either standard or essential supremum), with some weight x(t, x) (to be specified 
later), we are able to bound the terms containing u on both sides of the previous equation 
in such a way that 

IM| < 5(V)\\ X u\\ + C(f,g). (14) 

Then if S(V) < 1 we arrive at 

\M\<^^=-C{f,g,V). (15) 

In the case of pointwise bounds the \\x • || norm is a standard weighted supremum norm, 
which implies a pointwise decay estimate 

In the case of L°° norms, this inequality holds almost everywhere. 
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The condition 5(V) < 1 presents a bound on the strength of the potential \V\ and is 
a necessary assumption restricting the theorem. Although such a bound is not expected 
for non-negative potentials 1 the reason for the restriction of \V\ is that the technique of 
proof does not distinguish between negative and positive potentials. Therefore, everything 
what can be proved, must hold for negative potentials of the given strength \V\. Yet, 
negative potentials may have bound states (i.e. eigenfunctions of —A + V) which may grow 
exponentially in time, thus destroying the decay. Hence, the condition |V(a;)| < V / (x) k with 
k > 2 and Vq small enough must in fact avoid bound states. Indeed, this can be observed 
here, at least for spherically symmetric potentials where several bounds on the number of 
bound states are well-known p, Th. XIII. 9], e.g. the Bargmann's bound 



N ±l Wr)\*< * k (17) 



o 



or the Calogero's bound 



2 f°° 2 f\T 

N<-1 ^m d r<-^. (18) 



We get N = if N < 1, i.e. if V < (k - l)(k - 2) or V < ^(k - If. The bound, which 
appears in our proofs V < C~ k , although far from optimal, guarantees these conditions. 

We make an effort to estimate the numerical value of the constant C v ^ since in some 
applications the potential appears with a given value of Vq which is not (arbitrarily) small, 
e.g. waves propagating in the Schwarzschild geometry give rise to an effective potential. 
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II. MAIN RESULTS 
Definitions 

We denote by M + := (0, oo) and its closure by IR + := [0, oo). Analogously, we denote by 

M+ +3 := R + x M 3 and by 1+ +3 := 1+ x IR 3 . 

We denote by S(x,t) a sphere in M 3 with the center x and radius t, i.e. S(x,t) := {y G 
IR 3 : \y — x\ = £}, and by K(x, t) a (past-)cone in R^ +3 with the center at point (t,x), i.e. 
K(x,t) := {(s,y) E M+ +3 : \y — x\ = t — s}. The corresponding Lebesgue measures on S(x,t) 



There still does not exist an analogous decay theorem, there is only strong numerical evidence. In 
some weaker decay is proved for arbitrarily strong positive potentials. 
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and K(x, t) we denote by da(y) with y G S(x, t) and dn(s, y) with (s, y) G M x ]R 3 fi i), 
respectively. 

We introduce a frequently used short hand 

(x):=l + |s|. (19) 

We remind that the space L°°(M ri ) consists of Lebesgue measurable functions on W 1 having 
finite L°° norm which is defined as 

\\h\\ L °° := ess-sup \h(x)\ = inf{M : \h(x)\ < M a.e. on W 1 }. (20) 

Here a.e. on M n means almost everywhere on M. n in the sense that it does not hold at most 
on a set of measure zero in M. n with the standard Lebesgue measure on W 1 . For bounded 
functions h(x) and H(t, x) we will use the notation 

||/i||oo : = sup \h(x)\ and Halloo : — SU P \H(t,x)\ (21) 

x <& 3 (i,x)Gl+ +3 

identical in both cases if it does not lead to confusion (otherwise we will use the symbol sup 
explicitly). We define weighted-L°° spaces for measurable functions defined on space (IR 3 ) 

or space-time (M+ +3 ) 

heL™ & \\h\\ L «> := ||(-)^|U-(K3) < oo, (22) 



H G L~ & \\H\\Lft := \\(t+ \x\) r (t - \x\) p - r H(t,x)\\ LriW+)L „ m < oo. (23) 

The content of the last definition is that H(t,x) decays like l/t r along the lightcone and 
like l/t p or l/|x| p at timelike or spatial infinity, respectively. We will mainly consider spaces 
L™ p . We define some constants used throughout the article 

c p :=l/2(p-2), (24) 

°P := max {w^Yy 4 ) ' c v ] := max (ft 5 ) ' (25) 

C m := max(C«, C?) = max ( 9 _ , 5) , (26) 
C M := ^^^y max(2/(p - 1), 3). (27) 



Theorems 

First, existence and uniqueness of weak solutions is proved together with conservation 
of the energy. The theorem provides also L 2 -bounds on Vu(t) and d t u(t) which are not 
contained in the following decay theorems (because the derivatives there are considered only 
to exist in the distributional sense) and implies that in fact the first derivatives exist in the 
weak sense. 
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Theorem 1 (Existence and uniqueness) . With the assumptions of either Theorem Ha| or 
l2U the (weak) solution u of the wave equation (f2P-|iP is unique, belongs to C°(R+, H 1 ) fl 
C 1 (M + ,L 2 ) and gives rise to a conserved energy 

E[u(t)} := (V«(t), Vu(t)) + (u(t), Vu{t)) + (d t u(t), d t u{t)) 
= E[u(0)] = {Vf,Vf) + (f,Vf) + (g,g). [ j 

Moreover, Vu(t), d t u{t) £ L 2 for all t £ K+ u>/ia£ implies that the first derivatives exist in a 
weak sense, i.e. Vw(t), d t u{t) £ L] oc for all t £ R + . 

The case when additionally V > is somewhat simpler to treat and Thoe [3] has given a 
complete spectral theory including explicit representation of the evolution operator in action 
on the generalized eigenfunctions. It provides conservation of energy, but it does not seem 
to be of a big advantage when proving decay estimates of the type we are interested in. 

It is important to note that for continuous potentials there is no need of using such 
sophisticated tools like functional analysis for proving uniqueness of classical solutions. In 
this case the standard energy inequality is sufficient. Indeed, note first that V £ C 2 is enough 
to guarantee the regularity of classical solutions u £ C 2 and by differentiation of the free 
energy E (t) := ||Vw||^ 2 + IH|^ 2 with respect to time one obtains easily 

^ < 2V E (t), (29) 

what implies E (t) < E (0) exp(2Vot). This, together with the bound ||w(i)||x2 < ||u(0)||i2 + 
Jo ll^(^')IU 2 dt! 1 gives bounds on (u(t),u(t)) £ H 1 x L 2 in terms of the initial data. By the 
density argument these bounds remain true in case V, u £ C° (with weak derivatives and finite 
energy initial data), what guarantees uniqueness of solutions in C°(lR + ,i/ 1 ) fl C 1 (R + ,L 2 ). 
Unfortunately, this technique cannot be extended to the case of weighted L°° spaces, because 
C n are not dense in L°° . 

We state the results regarding the continuous pointwise and weighted-L°° estimates sep- 
arately. The first theorem contains a pointwise estimate for weak solutions which are only 
continuous and presents a generalization of the theorem of Strauss et al. jl[ for classical 
solutions (twice continuously differentiable). 

Theorem 2a (Decay in C°). Let (f,g) £ C 1 ^ 3 ) x C°(M 3 ) and V £ C°(M 3 ). If for some 
k > 2 

\V(x)\ < -% Va; £ M 3 with V < C: l k < 00 (30) 



(xy 



and for some m > 3 



l/(*)l<^pr> l v/ ^(^' l*^<Jp VxGR3 ^ 

then the wave equation (C]) with the initial data (TJ|) has a unique weak solution u £ C 
which satisfies 

W'-')i^ (t+W) ?- Wr . - v(^)^' +s < 32 > 



o/W 1+3 ^ 



with p := min(A;, m — 1) and 

1 — C P ,kVo 



n _ C m (fo + fi + 9o) /qq\ 
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The condition < fo/(x) m ~ l can be replaced by a weaker — > as \x\ — > oo. 

The pointwise estimate follows then from integration of |V/(x)| < f\/ (x) m from infinity to 
x along the radial direction. Therefore, it holds fo < Ji/(m2 m ). 

The second theorem generalizes the estimate further to weak solutions which belong only 
to a weighted-L 00 space. This seems to be the widest space compatible with the assumptions 
of spatial decay of the initial data and the potential. 

Theorem 2b (Decay in L™). If for some k > 2 

V := \\V\\ L - < C£ < oo, (34) 

and for some m > 3 

fo '■= ll/IUs-i < 00 h := ll V /IU^ < 00 9o ■= \\9\\l% < oo (35) 
then the wave equation (QP with the initial data |||) has a unique weak solution which satisfies 

|MU- < C, (36) 

with p := min(/c, m — 1) and 

n _ C m (fo + fi + go) , . 

i r- - !/ ' ^ ^ 

Note that the estimates of Theorem [2a] or [2b] 

|n(t,2)|<- ——, — — r <-V, a.e. for x G M 3 , Vi G M+ (38) 

1 v n ~ (t+ \x\)(t- \x\y- 1 ~ (x)p + \ J 

for p > 2 and some C(t) < oo for all t G IR + imply G L 2 for all £ > 0. 

Then, in a straightforward way we can prove the following, useful for further applications, 

Corollary 1 (Decay in presence of a source). Let the initial data satisfy 

fo ■= ll/IU-.i < oo, h ■= ||V/||i» < oo, g := \\g\\L°° < oo (39) 

with m > 3, the potential satisfy 

Vo := \\V\\l~ < C~l < oo (40) 
tuit/i k > 2 and the source F G L^° r satisfy for some q > 2 and 1 < r < q 

F := IK^FH^. < oo. (41) 
Tnen there exists a weak solution u(t, x) of the wave equation with potential and source terms 
Du + Vu = F, u(0,x) = f(x), d t u(0,x) = g{x) (42) 
and initial data u(0,x) = f(x), dtu(0,x) = g(x) which satisfies 

n n . Cm(/o + fi + 9o) + C rq F 

IMU « " T^v a • (43 » 

for p := min(/c, m — 1, r). 

If additionally (f,g) G C^M 3 ) nC°(M 3 ) ; V G C°(K 3 ) and F G C°(I+ +3 ) £nen n G C°(I+ +3 ) 
and £ne same estimates hold pointwise 2 . 



i.e. everywhere, instead of almost everywhere. 
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Corollary 2 (Decay of derivatives). The weak solution of the wave equation (TJP with the 
initial data |I)) satisfies 

u, \Vu\, \V 2 u\,d t u, d 2 t u G L™ p (44) 
with p := min(/c, m — 1) provided 

f, |V/|, |V 2 /I, |V 3 /I, |V</|, |V 2 ^| eC m > 3 (45) 
V, | W|, I W| G fc > 2 and \\V\\ L f < C~l < oo. (46) 

If additionally {f , g) G C 3 (M 3 ) nC 2 (R 3 ) and V" G C 2 (M 3 ) fenu6C 2 (R + ) and the estimates 
hold pointwise 2 . 

Here, we have introduced notation |V n /i| := Ylai a»=i \9 ai ---d an h\. This is the same result 
(in the u G C 2 case) which was obtained by Strauss and Tsutaya in [l|. 



III. BASIC ESTIMATES 



The first lemma estimates an integral appearing in the Duhamel representation formula 
containing initial data. It was proved in somewhat different form by many authors [l|, @, 0] , 
but here we prove only a simplified version, which is needed for our goals. Actually, this 
estimate can be still improved by one power of (t — \x\), what better applies to long range 
initial data, which are not of our main interest. 



Lemma 1. For t > 0, p > 2 and 



"=±t ^ (47) 



^Js(x,t) (y) 



it holds 



Moreover, for p > 2 



and for p > 3 



l < ' > , (;- ,)^- (48) 

j_ r < ^ (49) 



My) < cf 

47Tt 2 J S (x,t) {VY' 1 ~ (t+\x\)(t-\x\)P-*' 1 ] 

The next lemma estimates an integral appearing in the representation formula containing 
a source term. We follow the proof of Georgiev et al. Similar calculation has been done 
by Strauss et al. in [lj and earlier by John [9|. 

Lemma 2. For t > 0, x G M 3 and 

j. = ]_[ dn(s,y) 

4vr J K{x>t) (t - s)(y)o(s + \y\)(s - M)*- 1 



for q > 2 and q > p > 1 it holds 



1/1 < 9lA . (52) 

1 1 " (t+ \x\)(t- \x\y- 1 V } 
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Next, we prove the general representation formula for distributional and weak solutions. 

:mma 

formula 



l-|-3 

Lemma 3. For any distributions f,g G P'(1R 3 ) and F G V(R 1+3 ) supported on R, the 



S(t- \x-y\) 3 f 5(t-\x-y\) 



«(«. *) = ft / ' 4rt "W + / 4 ; f "' <ta)*v 



3 5(t-S- 



(53) 



1 _l_3 

defines a distribution v G 'D / (IR 1+3 ) supported on IR + which satisfies the wave equation 

Dv = F, v(0,x) = f(x), d t v(0,x) = g{x) (54) 
in the weak sense, i.e. for any test function tp G C^(M 1+3 ) 



dt I d 6 x Dtp(t,x) v(t,x) = - I d A x d t ip(0,x) f(x) + / d A x ip(0,x) g(x) 

J J (55) 

+ dt d 3 x tp(t,x) F(t,x). 



If f,g,F are ordinary functions (i.e. regular distributions) or in L\ oc and v G L\ oc is defined 
almost everywhere in R + by the Duhamel's formula 

(56) 



4?r Js(x,t) t 4vr Js{X)t) 

+ ir ^yl dK(S:V) , 



4vr JK(x,t) t - s 

then v is a weak solution, in the sense specified above, of the wave equation (54\ ). 

The following lemma, an estimate for solutions to the free wave equation with prescribed 
initial data, has been proved by Strauss et al. [l[ for classical solutions, i.e. for (/, g) G 

C 3 (IR 3 ) x C 2 (M 3 ) leading to u G C 2 (l+ +3 ). We weaken the assumptions to (f,g) G C^IR 3 ) x 

C°(M 3 ) and obtain u G C°(I+ +3 ). 

Lemma 4a. Let the data (f,g) G C X (1R 3 ) x C°(1R 3 ) and satisfy 

l/(*)l<^pi. W^^p |y(x)l -<Jp VxGR3 (57) 

for some m > 3. Then there exists a unique weak solution v(t, x) = Io(f, g) of the free wave 
equation 

Dv = 0, v(0,x) = f(x), d t v(0,x) = g(x). (58) 

X_l_3 

Moreover, it is continuous in [t, x) G M + and satisfies 

\v{t,x)\<- , C y ,d ] ,. , V(t,x) Gl+ +3 , (59) 

1 v n ~ (t+ \x\)(t- \x\) m ~ 2 y ' + ( ; 

where C(f, g) := C m ■ (g + f\ + / ). 
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As next, we further weaken the assumptions to V/, g G and / G -k^-D i- e - to 
weighted L°°(IR 3 ) spaces. 

Lemma 4b. Let the data (f, g) G L'^_ 1 x with m > 3 satisfy 

fo ■= Wfh^ < oo, /i := ||V/||l» < oo, g := ||#|| L - < oo. (60) 
Then there exists a unique weak solution v(t,x) = I (f,g) of the free wave equation 

Dv = 0, v(0,x) = f(x), dtv(0,x)=g(x) (61) 

which satisfies 

IMU^ < C(f,g) := C m ■ (g + fi + /„). (62) 

The following lemma, an estimate for solutions to the wave equation with source, has been 
proved by Strauss et al. [l| and Asakura [s| for classical solutions, i.e. for F G C 2 (lRi +3 ) 



-y , i ,g 1+3 

leading to u G C 2 (R + ). We weaken the assumption to F G C°(1R + ) and obtain u G 
C°(I+ +3 ). 

1 _l_3 

Lemma 5a. Let the source F G C°(lR + ) and satisfy for some q > 2 and 1 < p < q 

\F(t, x) | < , ... F °, ,. . V(t, x) G I+ +3 . (63) 

T/ien i/iere exists a weak solution v(t,x) = Lq(F) of the free wave equation with source 

Dv = F, (64) 

1+3 

and null initial data v(0,x) = 0, d t v(0,x) = 0. Moreover, it is continuous in (t,x) G M + 
and satisfies 

I'ft'^ ^i^V- v ^ )eS + + - (65) 

1+3. 

As next, we weaken the assumption to F G (x)~ q L^° p , i.e. a weighted L°°(IR + ) space. 
Lemma 5b. Let the source F satisfy for some q > 2 and 1 < p < q 

F := \\(x)*F\\ L ~ < oo. (66) 
T/ien there exists a weak solution v(t,x) = Lq(F) of the free wave equation with source 

Dv = F, (67) 
and null initial data v(0,x) = 0, dtv(0,x) = 0. Moreover, it satisfies 

IMIl- < C M F . (68) 



10 



IV. EXISTENCE AND UNIQUENESS IN THE ENERGY SPACE 



For the free wave equation with V = it is well known that the homogeneous Sobolev 
spaces H S (M. 3 ) x H S ^ 1 (M 3 ) have the property that the evolution operator W(t) generating 
the solution (u(t) , d t u{t)) = W(t)(f,g) of the free wave equation is unitary and the norm 
defines a naturally conserved energy 



U W% S + ll<9t«ll!- s -i 



E s [u(t)} : = 
with the case s = 1 of our special interest 

E Q [u(t)} := J (|Vu| 2 + |«9H 2 ) d 



x. 



(69) 



(70) 



This is no more true in presence of the potential V . Yet, if the potential satisfies some 
weakness conditions we show that there exist perturbed homogeneous Sobolev spaces 
7^ S (]R 3 ) x 7^ S_1 (IR 3 ) in which the evolution is unitary 

We start from the observation that initial data (/, g) satisfying the decay bounds (|3T|) or 
(|35|) of Theorem [2a] or [2b] respectively, 



\f{x)\ < 



/o 



(x) 



m-1 ' 



|V/(x)| < 



fi 



(x)> 



\g(x)\ < 



go 

(x) n 



m > 3, 



(71) 



holding either everywhere or almost everywhere in M 3 , belong to the following spaces: Xo : = 
H 1 x 1? (the free energy space), X 1 := H 1 x L 2 and X E := Ti s x 1-L S ~ X (defined below) with 
s — 1, which we will call the (true) energy space for the problem at hand. In the following 
we will prove that if the potential satisfies the bound f[3"Uj) or of Theorem [2a] or l2bl 
respectively, 

<t|^ (72) 

again everywhere or almost everywhere in M 3 , then there always exists a unique solution 
(u, u) G C(R+, X ), C(M+, X\) or C{R + ,X E ) of the wave equation ([!]) for any initial data 
(/,<?) £ Xo,Xi or Xe, respectively, but only in the true energy space X E the evolution is 
unitary and the corresponding energy conserved when V ^ 0. 
In dimension n = 3 the bounds (171]) and (172]) imply 



Moreover, from (172]) with k > 2 follows 



and 



In 



v e L . 



< 



\x\ 



(l + |x|) 2 

and with (•,•), being the standard scalar product in L 2 , we have for any h £ i/ 1 

|Mx)| 2 



\(h,Vh)\ 



V(x)\h(x)\ 2 d'x 



< Vr 



J l ^~ d3x ^ 4Vo J \^h\ 2 d 3 x = W \\Vh\\ 2 L2 



by Hardy's inequality for n = 3. For any h £ L 2 we have 



1(^^)1 



y(x)|/i(ar)|Vx 



< iiy| 



l^lli- < Vo 



2 

L 1 ■ 



(73) 



(74) 



(75) 



(76) 
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For the classical solutions u G C 2 



we know that the energy, defined by 



E[u(t)\ := J (\Vu\ 2 + \d t u\ 2 + V\u\ 2 ) rf 3 x, 



(77) 



is conserved in evolution, what can be proved by differentiation with respect to time. This 
technique does not work for weak solutions, for which the derivatives exist only as distribu- 
tions. For our rough data (/, g) we have from (1751) 



\E\u 



J (\vf\ 2 + \g\ 2 + v\f\ 2 )d' 



x 



<l|V/||i 2 + |b||! 2 + y ||/|li 2 <oo. (78) 



This energy is positive definite for potentials which are weak enough 3 and satisfy Vq < |, 
what follows from (1751) and 

\\Vf\\h + <f,Vf)> ||V/||i a -4V j \Vf\ 2 d*x = (l-Wo)\\ V/Hi, > 0. (79) 

In this section we will show that there is a unique evolution in the space of solutions with 
finite energy. 

First, consider A := — A + V defined on C^°(IR 3 ). For V G L°° H L 2 , what is guaranteed 
by (172|) . the operator A is essentially self-adjoint on C^°(1R 3 ) and has a unique self-adjoint 
extension (which we denote for simplicity again by A) to T>(A) =T>{— A) = {/ G L 2 : Af G 
L 2 } = H 2 (see e.g. 0, Th. X.15]). 

For / G H 2 we can partially integrate and find as in (179"1) that 



(/, Af) = (/, -Af) + (/, Vf) = || V/||i 2 + (/, Vf) > (1 - 4Vb)j|V/||| 2 > 0. 



(80) 



Moreover, we have, again by Hardy's inequality, 



(/^/)=nv/iii 2 +a,F/)>iy 



i/(^)i 



\x 



> 



4 b 



In 



U1 



^ ^ > 5 



rf 3 x + y ik(x)| i/(x)i 2 d 3 x 
I/WI 2 



(81) 



drx 



\x\ 



where S :— \ 



V > 0. Then, (/, Af) = implies / ^f - d 3 x = andhence \ f{x)\ = 



(a.e. on IR 3 ). This shows that A is positive definite. Then, according to [Uj, A has a unique 
positive definite self-adjoint square root B := A 1 ! 2 with T>(B) = {/ G L 2 : Bf G L 2 } D 
P(y4). In spirit of the spectral theorem we can define positive self-adjoint A s l 2 for any s > 0. 

We recall that the homogeneous if s and inhomogeneous H s Sobolev spaces are defined 
as completions of w.r.t the norms 



\h\ 



-Ay/ 2 h\ 



L 1 



and 



\h\\ 2 
\ u \\h s 



ElNA* = ElK- A )* /a *. ( 82 ) 



fc=0 



fc=0 



respectively, and have the property H C\H S = H s . In order to use the standard construction 
of a unitary evolution given by the self-adjoint generator A we define perturbed homogeneous 



It is also positive for all V(x) > 0, but here we are interested only in bounds of the form ([7^)1 



12 



H. s and inhomogeneous 7i s Sobolev spaces as completions of w.r.t the norms (in case A 
is positive but not positive definite || • ||^ s are only seminorms) 

s s 

\\h\\ ns := \\A s ' 2 h\\ L , and \\hf w := £ \\h\\\ k = £ \\A k ' 2 h\\ 2 L ,. (83) 

fc=0 fc=0 

for sgR In our situation A is positive definite, so || • ||^ s and || ■ ||^ s are norms and lri s 
and TC S are (complete) Banach spaces. 

Now we relate Ti, 1 to energy. Consider any / e C^°(IR 3 ) for that we can partially integrate 



(A 1 / 2 f,A 1 / 2 f) = (f,Af) 

(f, -A/) + (/, Vf) = (Vf, V/) + (/, Vf) =: ll/H Jj, 



(84) 



The last expression defines what we will call "partial energy" norm (its norm properties 
are inherited from || • H^i). If we define the "partial energy" space Ye by completion of 
C^°(1R 3 ) w.r.t. this norm then the spaces Ye and H} are identical, because both are defined 
as completions of C|?°(R 3 ) w.r.t. to the norms which coincide 4 on the whole C^°(IR 3 ). Hence, 
Ye = Ti 1 (see also 12j for similar considerations). 

We can now define the energy space 5 as X E '■= Y E x L 2 = H 1 x L 2 or equivalently as a 
completion of C^°(1R 3 ) w.r.t. the norm 

Uf,9)\\ 2 x a ■= llV/lli, + (f,Vf) + \\g\\ 2 L2 = \\f\\% + \\g\\ 2 L2 . (85) 

Xe is a Hilbert space with the scalar product 

((fu9i)\(f2,92))x E := (A 1/2 fx, A 1/2 f 2 ) + ( 9l ,g 2 ), {fugi) 6 H 1 X L 2 , (86) 

what follows from the fact that every Banach space H, s can be made to a Hilbert space 
equipped with the scalar product 

(;-) H8 :=(A s / 2 ;A s / 2 .) L2 . (87) 

We enhance the facts that A and B are self-adjoint on L 2 with the domains T>(A) = 7i 2 C L 2 
and T>(B) = Ti 1 C L 2 , respectively, by the observation that i? defined on H 1 is self-adjoint 

with Vfr(B) = n 1 nn 2 c n 1 . 

For more flexibility in the applications we also relate the homogeneous and inhomogeneous 
Sobolev spaces. For s = we trivially have H° = H°. From fl75]) and $72} it follows 

< (1 - 4T/ )||V/i|| 2 2 < HV/illl, + (/!, ^) < (1 + W )\\Vh\\ 2 L2 , (88) 



4 More precisely, the completions are defined by Cauchy sequences in C5°(K 3 ) w.r.t. the given norms. Since 
both norms arc equal on (^(M 3 ), the Cauchy sequences w.r.t. to the two norms are identical. 

5 We choose Cartesian products instead of direct sums in order to avoid confusion with norms. A direct 
sum of Banach spaces induces a norm and a direct sum of Hilbert spaces induces an inner product, which 
in turn induces a norm, but these two norms are different. Since we treat the same spaces as Banach 
and Hilbert we use rather Cartesian products and specify the norms explicitly every time we define a new 
space. 
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what means that there exist constants C\,Ci > such that < ||/i||y^ < C^H^Hjji, 

hence the norms in H 1 and Ye are equivalent and H 1 = Ye as sets. Since Ye = Ti. 1 we 
have H 1 = 1i} as sets with equivalent norms. The same is also true for the inhomogeneous 
spaces, namely H 1 = Ti 1 as sets with equivalent norms. 

Furthermore, since || V/i||^2 < Vo||/i||l2 we have on the one hand that h G L 2 fl H 2 implies 
\\Ah\\ L 2 < || - Ah\\ L 2 + Vo||/i||L2 < oo and hence h G H 2 . On the other hand, h G L 2 HH 2 
implies || - Ah\\ L 2 < \\Ah\\ L 2 + Vo||/i||i,2 < oo and hence h G H 2 . That gives H 2 = H 2 . 

Now, we are ready to formulate and prove global existence, uniqueness and energy con- 
servation for solutions of the wave equation in the energy space Xe- 



Theorem 1'. Assume the potential V satisfies the bound ( f7l?| ) with Vq < 4. 



a) For all initial data (f, g) G Ti. 1 xL 2 there exists a unique weak solution u G C°(1R + , li, 1 )^ 
C 1 (1R + ,L 2 ) of the wave equation (Qp. Moreover, the energy defined by the norm in Xe is 
constant in time 



\\(u(t),d t u(t))\\x E = (V«(t), Vu(t)) + (u(t),Vu{t)) + {d t u(t),d t u(t)) 
= (Vf,Vf) + (f,Vf) + {g,g). 



(89) 



b) Initial data (f, g) G H 1 x L? or H 1 x L? give rise to a unique solution u G C°(R + , H 1 ) fl 
C 1 (1R + ,L 2 ) or C°(W + , H 1 ) nC 1 (M + ,L 2 ), respectively, but the corresponding norms are not 
conserved in evolution. 

Proof, a) In order to prove existence of a unique (weak) solution in Xe = 'H 1 x L 2 we take 
advantage of the functional analysis in Hilbert spaces. We show that the evolution in Xe 
is generated by a self-adjoint operator M and use the famous theorem of Stone to show 
that there exists a unique strongly continuous group of unitary operators W(t) defining the 
evolution of u. Unitarity of W(t) guarantees conservation of the norm in X, which represents 
the energy of u(t). 

First, we write the wave equation (TjQ) as 

d 2 

- —u = Au = (-A + V{x))u. (90) 

A with T>(A) = Ti 2 is self-adjoint on the Hilbert space Z := L 2 and positive definite when 
the bound ^ is satisfied with V Q < \ (cf. ((79])). Its unique positive definite square root 
B := A 1/2 is self-adjoint with V Z (B) = H 1 C Z on the same Hilbert space Z = L . On the 
Hilbert space Y := Iri} with the scalar product (•, -)y := (B-, B-)z the operator B is again 
self-adjoint with V Y {B) =H 1 nH 2 CY. 

Next, the wave equation (|9"U|) can be written as a system of first order differential 
equations in time on X Z z '■— T>z(B) x Z = H 1 x L 2 with ((/i, gi), (/ 2 , g2))x zz ■= 
(fx, f 2 )z + (Bf x , Bf 2 ) z + ( 9l , g 2 ) z 

^-U(t) = -iMU(t), U(0) = ( f ) (91) 



dt w w y ' \g 

with 

U{t) := ( ^ \eX zz and M := z ( _ ° J). (92) 
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M generates a strongly continuous semigroup Wzzif) on Xzz which solves the above ini- 
tial value problem, but M is not self-adjoint on Xzz (not even symmetric) and hence 
Wzz(t) is not unitary on Xzz- Therefore, we consider now a completion Xyz of Xzz in the 
norm \\(f u g\)\\ XYZ := ||-B/i||| + IM|| = ||/i||y + IM|| equipped with the scalar product 
((fi,9i),(f2,92))x YZ ■= (Bf 1 ,Bf 2 } z + {9i,92)z = {h,h)y + (9i,92)z- It turns out that 
X YZ = Y x Z = H 1 x L 2 = X e . Then the closure of M with V(M) = V(A) x V Z {B) = 
7i 2 x 7i l is self-adjoint on Xe (see 13, Lem.7.7] for a general proof or [3] for the case V > 
treated more explicitly) and generates a strongly continuous unitary group W(t) on Xe 
given by the unique extension of 

to Xe (see also [h], sec. X.13] or 14, sec. 6.2] for similar considerations). Such W(t) defines 
a strongly continuous evolution in X E 

m -wwV>-wM(' a ). (94, 

Thus, we have U = (u,u) G C°(M + , H 1 x L 2 ) or, in other words, u G C°{W + ,H 1 )nC 1 {W + , L 2 ). 
From the unitarity of W(t) we get — ||^(0)||x B what expressed in u(t) reads 



2 i II - /+MI2 || t\\2 , || 1 1 2 



U (t + U 



Til 1" || "Will, 2 - IIJ llTi 



Mil- (95) 



Since the above norm in H} x L 2 is equal to the energy norm || ■ \\ Xe we get the energy 
conservation fl89|) . 

b) Conservation of energy in Xe implies that the solution also belongs to the free energy 
space X = H 1 x L 2 

E [u(t)] = ||Vu(t)|| 2 2 + \\d t u{t)\\ 2 L2 < dE[u{t)] = dE[u(0)] < oo (96) 

with d := (1 - Wo)' 1 > 1 what follows from ([79]) and V < \ (for V > the same holds 
with C\ = 1). Moreover, by (1751) . 

E [u(t)} < C x E[u{G)\ < C 2 E [u(0)}. (97) 

Summarizing this, there is a continuous (but not unitary) evolution (/, g) G H 1 x L 2 — ► 
(u(t), <%«(*)) G if 1 x L 2 for all t > 0. 

In case when additionally / G L 2 we get u(t) G L 2 for all t G R + , because 



l«(*)ll 



L 2 



m(0)+ / d t u(t')dt' 



< \\f\\ L 2 + / ||5 t «(t')|U^t' 



(98) 

<\\f\\ L 2+ [ ^E[u{t)]dt' < \\f\\ L 2 + v^KO)] t < oo Vt>0, 



and a continuous evolution (/, g) G if 1 x L 2 — > d t u(t)) e H 1 x L 2 for all t > 0. □ 
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Theorem [T] 



Proof. In ()73j) - (j78p we have shown that the assumptions of Theorem [2a] or [2b] imply finiteness 
of the energy of initial data i?[«(0)] < oo. Since by f[85| 



\g\\ 2 L2 = E[u(0)} < oo, (99) 



we have (/, g) G H 1 x L? . The condition V§ <C \ implies Vb < \ for any p > 1, k > 2. So 

we can apply Theorem [Fl a) and obtain a unique solution u G C (M + , H 1 ) PI C 1 (M + , L 2 ) with 
conserved energy 

E[u{t)] = \\(u(t),d t u(t))\\x E = (V/, V/) + (f,Vf) + (g,g) = E[u(0)}. (100) 

Since the assumptions of Theorem [2a] or [2b] imply also (/, g) G if 1 x L 2 we have by Theorem 
QZ]b) that u G ^(1+,^) nC 1 (I+,L 2 ). 

This gives_u(t),Vu(t),d t u(t) G L 2 for all t e R + what implies u(t),Vu(t),d t u(t) G Zq 00 
for all t G R+, i.e. existence of the first derivatives in the weak sense, because for any 
compact set Q C M 3 and h = u(t), Vw(t) or dtu(t) Holder's inequality gives 

Jjh\dfi< (J^dfj) ' -(Jjh\ 2 d^ ' < /x(^) 1 / 2 • < oo. (101) 

The last step is to show that such constructed u G C°(R + , H 1 ) n C 1 (M + ,L 2 ) really solves 
the equation ([T]). Indeed, it solves the equation in the distributional sense (see [HI, Th. 3.2] 
for a more detailed discussion). Since additionally u G L\ oc it solves §S§ in the weak sense 

a. □ 



V. PROOFS OF THE ESTIMATES 
Lemma [TJ 

Proof. For t > we have 

" 4vr ./^t) 4vr 7 S(0)1) + tu;|) 4vr 7 S(0)1) + iu;|) 

Using the fact that \ou\ — 1 and introducing the angle 9 between x and u we get 

\x + tuo\ 2 = \x\ 2 + t 2 \cu\ 2 + 2t(xu) = \x\ 2 + t 2 + 2i|z| cos6> (103) 
hence in polar coordinates (9, <fi) 



,2 rir f2ir 

I = — d9sin9 ' 



4vr J J f(^/\x\ 2 + t 2 + 2t\x\cos9) 

t 2 r +1 J 1 



(104) 



2 i-i /(Vk| 2 + t 2 + 2t|x|x) : 



1(3 



where we substituted x '■= cos 9. Now, changing variables A 2 := \x\ 2 + t 2 + 2£|x|x and 
\d\ = t\x\dx we get 



_ t 2 f t+ ^ XdX 



t+\x\ 



XdX 



2 J lt _ {xll t\x\f(X) 2\x\J ]t _ M f(XY 
Next, choosing /(A) = = (1 + A) p and p > 2 we have 



AdA 



2|x|,/ |t _ w| (1 + A)f" 
This simple integral can be evaluated exactly for p — 2 > and then estimated 



1 



1 



2\x\ (p-2) (1 + A)f- 2 
1 



p- 2 



(p -1)(1 + A) 



|t+|x|| 



\t-\x\\ 



t 



+ 



2|x| L(p-2)(1+ \t- \x\\)p~ 2 (p-2)(l + \t+\x\\)p- 2 
1 1 



(p - i)(i + |t + i^idp- 1 (p - 1)(1 + \t- \x\dp- 1 _ 



<o 



< 



2(p-2)]x| (t- |x|)P" 2 ' 



It gives the estimate ( l48|) with c p := l/2(p — 2). 

We proceed further in two ways: for q > 2, on the one hand we have from (11071) 



1 



S(x,t) 



d<r{y) 

(y) q 



< 



\x\(t- \x\)i~ 2 ' 



and the other hand, taking the supremum out of the integral, we find 



1 

4vrt 



S(x,t) 



dcrjy) < J_ 
{y) q ~ Ant 



sup — 

\ y -x\=t (y) q 



1 Airt 2 _ t 

W) a[y) ~Art{t-\x\y ~ {t-\ x \y 



Combining both inequalities H108[) , (I109[) for q = p > 2 we get 



1 

4vrt 



da(y) 



< 



> S (x,t) (y) p ~ (t-\x\)p- 

By little algebra it can be shown that 



mm 



2 Cn 



2\ x y (t- \ x \y 



2 Cr, 



< 



9 c Ti 



2\x\ ~ (t + 

t 4 

< 



(t- |x|) 2 " (t+ \x\) 



< 



(t-\x\) 2 - (t+\x\) 



for \x\ > -, \x\ > - 
4 4 

for Ixl < -, t > 1 
- 4' - 



for \x\ < -, t < 1. 

1 1 - 4' - 



(105) 



(106) 



(107) 



(10£ 



(109) 



(110) 

(111) 
(112) 

(113) 
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Therefore, for p > 2 



where 



da(y) 



< 



a 



(i) 



*rtJs(x,t) (y) p ~ (t+\x\)(t- x)p- 



Cjp := max ( 9c p ,4, 



9 



max 



9 



4 J \2{p-2Y 
Now, combining inequalities (llQ8j) . ( 1109[) for g = p — 1 > 2 we get 



1 



4vrt 2 



< 



mm 



S(x 



t) (y) P -i- (t-\x\)r* \2\x\t> (t-\x\y 



2c„_i 



By some more algebra it can be shown that 



1 \x\ 

for - < — < 2 and \x , t > 1 : 

2 - t - II'- 

M M i 

for — > 2 or — < - : 
t ~ t ~ 2 



for \x\, t < 2 : 



So we have 



mm 



2c„_i 



< 



max(6 c p _i, 5, 5) 



2|z|t ' (t- |x|)V ~ (t- \x\)(t+ \x\)' 



So we finally get for p > 3 



do-(y) 



< 



(2) 



where 



C^ 2 ' ) := max(6cp_i, 5) = max f — - — ,5 J 

VP - 1 / 



(114) 



(115) 



(116) 





2 c p _i 


< 




6c p _i 


(117) 




2|x|t 


(t~ 


\x\){t+ \x\) : 




1 


< 




5 


(118) 


(* 


- M> 2 


(t~ 


\x\)(t+ \ x \y 




1 


< 




5 


(119) 


<* 


- \x\) 2 


(t~ 


\x\)(t + \x\) 



(120) 



(121) 

(122) 
□ 



Lemma [2] 

For completeness we cite the proof after Georgiev et al. 

Proof. For t > and the continuous integrand in / we use a version of Fubini's theorem and 
split the integration over the cone K(x, t) into integrations over time and over spheres 



i r l 

I=— ds 
4W 



da(y) 



1 



S{x,t-s) 



(t-8)(y)*{8 + \y\)(8-\y\)r-i- 



(123) 
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rt pt—s+\x\ 

ds I d\ 



For integration over the spheres S(x, t — s) we use lemma [T] and obtain 

A 

2M J ~" J\^ s - [x f" (\)i(s + \)(s-\)p-i 

1 ft ft~8+\x\ 1 

/ ds I d\ 

2N7o J\t-s-\ x \\ (l + A)^- 1 (l + s + A)(l + | S -A|)f- 1 



(124) 



One can estimate the integrand for q > p 
1 



< 



3 9-l 



(l + A)"-i(l + s + A)(l + |s-A|)P- 1 " (l + s + A)f V(l + I«- AD 9 " 1 (1 + A)^- 1 

(125) 

because for < s < 2A it is smaller than the first term and for s > 2A smaller than the 
second one (see jEf or [l| for more details). Now, changing variables a := X + s and (3 := A — s, 
we get for q > 2 



/ < 



3 q- 



1 /-i+M 



2 b 



l*-MI 



(1 + a) 



+00 



d(3 



1 



(1 



9-1 



|g±g| 



q-l 



< Vr" (1 + 2 1 -") / , , . 

" 2|x|g-2^ . J sJ Mx]l (l + a )p 

<3/2 



(126) 



The last integral can be estimated either by taking supremum of the integrand or by inte- 
grating it explicitly 

j_ / <+N da r (1 + |t - ixiD-p- 1 1 f (t - Ix))- 1 

• 2|x| 7| t _ N || (1 +a)P-\ (2(p - ^kD-Hl + I* - |x||)" p (t - | x |>p-i \ (2(p - l)!^)" 1 

(127) 

It can be easily shown that 



< 2 ^ p , li f or 2\x\ > t, 2\x\ > 1 



2p\x\ ' (t + \x\) 

1 3 

< 



(t-\x\) ~ (t+\x\) 

1 3 

< 



and hence 



Finally, we obtain 



with 



(t-\x\) ~ (t+\x\) 

max(2/(p-l),3) 



for 2|x| < t, t > 1 
for t, Ixl < 1 



J < 



J < 



G 



a 



(t+ \x\)(t- \xDp- 1 

3 



P.9 •" 



2 (<z - 2) 



max(2/(p- 1),3). 



(128) 
(129) 
(130) 

(131) 
(132) 
□ 



19 



Lemma [3] 



Proof. For any distributions f,g G D'(IR 3 ) and F G T>'(K 1+3 ) with support on IR^ 3 the 
representation formula ( l53l) can be written as 



47r(t - s) 

+ / dsi d * y ^- s -\*-y\ 



0, 



*(•*-!■» I) 



+ 



4-7T-/ 



47r(t - s) 
**(■*)/(■«) 



e(s)F(s,y) 

*(•*-!-. I) 



(t, x) + 



4?rv 



«(-t)»(- a 



4vr(t - s) 



(t,x) 



W(-t)f(- x ) + S(- t )g(- x ) + 9(- t )F} 



(t,x), 



i.e. a sum of M 1+3 -convolutions of distributions 6(t)f(x),6(t)g(x),6(t)F(t,x) G £> 



(133) 



l-|-3 

with <5(t — |x|)/(47rt) G T?'(IR 1+3 ), all supported on R + . The formula defines again a 



^1+3 



distribution u G P'(R 1+3 ) [16j supported on M^_' In order to show that v solves weakly the 
wave equation ( 1541) . i.e. v satisfies ( 1551) . we transform the l.h.s. of ( 1551) for any test function 
V? G Cg°(R 1+3 ) 



(it / d 3 xv(t,x) D(p(t,x) 



8{-t-\-*\) 

47T- t 

Sh-\- X \) 
47W 



(it / d 3 x 



*[S'(- t )f(- x ) + 5(- t )g(-*) + 0(-t)F] 



(t, x) Oip(t, x) 
(t,x) <p(t,x) 
(t,x) (p(t,x) 



(134) 



where we have first used the definition of differentiation of distributions (corresponding to 
partial integration) and then the theorem on differentiation of convolutions [l6[ ■ Now, with 
the well-known distributional identity 



we arrive at 

dt I d 3 xv(t,x) Oip(t,x) = 

dt J d 3 x [5'(t)f(x) + 5(t)g(x) + 9(t)F(t,x)] <p(t,x) 

d 3 x f(x) d t ^(0,x) + / d 3 x g(x) <p(0,x) + dt I d 3 x F(t,x) <p(t,x) 



(135) 



(136) 
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In case when f,g and F are ordinary functions or in L\ oc , then they are also (regular) 
distributions with, in general, distributional derivatives. For them the statements above 

remain untouched. If additionally v G L\ oc is defined almost everywhere as a function 

1-L3 

on M + by the Duhamel's formula (156|) . which is equivalent to (1531 with only alternative 
notation for the measures (5(...)d 3 y — > do~(y)), then it also solves weakly the same wave 
equation, i.e. its derivatives exist in the distributional sense and v G L\ oc satisfies (j!36p . □ 



Lemma I4al 

Proof. Let v(t, x) for t > 0, x G M 3 be given by the formula 

v(t,x) = — ——da(y) + — da(y) 

^.IS(x,t) t 47T J S ( x ,t) t 2 



— g(x + tuj) da(uj) + — / u ■ Vf(x + too) da(uj) (137) 

47r Js(0A) 4tt Js(oa) 



'S(0,1) ^" ^5(0,1) 

f TZ I f(x + tuj) da(u). 

'5(0,1) 



47T 



Since g, /, V/ G C° these integrals exist for all (t, x) G M^ +3 and are finite. Moreover, since 
all three above integrands are continuous functions of (t, x) and are uniformly bounded by 
9o, fi, fo £ -^ 1 (5'(0, 1)) for all (t,x), from continuity of Lebesgue integrals, we obtain that v 
is also continuous in (t,x) G M 1 ^ 3 . The first two integrals are bounded by 47rg and 47r/ 1; 

respectively, hence limt^ + 2) = /(#) = f(0,x), so we have v G C°(1R + ). Then also 
v G L\ oc and from lemma [3] it follows that w solves weakly the wave equation fl58l) . 
Therefore, we are able to apply the pointwise estimates 



g(y)\ , |V/(y)| , 



t t t 



2 



da(y) 



< ll(-> m g||oo + IK-) m V/llco / da(y) IK-r-Vlloo f da(y) 

J S(Xtt) (y)»» + 4^ ./ 5(a||t) ^ ^ 

S'o + Zi /" da(y) f f da(y) 



with 



4vrt J s{x ,t) (y) m 4vrt 2 (y)™" 1 
From lemma [U eq. (1491) and (|50|) . for m > 3, we finally obtain 

l«^fflFB^ v(i ' l)eR ' +3 (139) 

C„ := max(C«, C<?) = max (" ^J* 5) . (140) 

The inclusion of the boundary t = is trivial, because from the bounds on the initial data 
we have 

|i>(0,s)| = < ^#-r < ° m ; ^ f \ + /°\ VxGM 3 . (141) 

1 v Jl lJ v ;| _ (M)™- 1 ~ (0 + \x\)(0- \x\) m ~ 2 K ' 

□ 
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Lemma I4bl 



Proof. We start again, as in proof of lemma Hal, with the formula 

v{t, x) = t~t [ 9(V) Mv) + / [(y - x) ■ Vf(y) + f(y)} da{y) (142) 
4vrt Js( x ,t) 4irt z J s( t) 



for t > 0, x G M 3 , but now we cannot use the pointwise estimates on / and g. We must show 
that the integrals over the spheres make sense almost everywhere and that the resulting v 
belongs to Lf m _ v 

For the moment let's forget the function / and concentrate on the first integral only. 
We can cover R 3 with spheres S(x,t) of same origin x and various radii t G M + , i.e. 

3 , which itself plays no role in integration and 
\\(-) m g\\L~ = go < oo we have \g(x)\ < go/(x) m 
a.e. on R 3 . The set 

A:={x:\g{x)\>g /(x) m } (143) 

has therefore measure zero in R 3 . From a variation of Fubini's theorem 6 for sets it follows 
that for almost all t G R + the measure (on the spheres S(x,t)) of A\s( x ,t) is zero. In other 
words, for almost all t G R+. and all x G M 3 it holds 



]R 3 ^ ^(x, •) x R + up to one point x G 
measurability. Since g G L^ 1 and ||(?||l^> 



k(y)l< Wd/>' 



(144) 



for almost all ?/ G S(x,t), what can be written as || (■) m g\\L ao (m + )L°°(S{x,t)) < fi'o- 

Since jel™ then also G L°°{R 3 ) and G L 1 1 oc (M 3 ). t From R 3 = S(x, •) x R + and a 
variation of Fubini's theorem 6 it follows that for almost all t G R + and all x G R 3 we have 
flisfot) £ -^ 1 (5'(x, t)) (we skip "local" because S(x,t) is compact inM 3 ), i.e. |0(y)| ^cr(y) 

exists and is finite. Therefore, for these (t,x), using (j!44p . we can estimate 

da(y) 



g(y) da(y) 



S(x,t) 



< 



\g(y)\ da(y) < g 



S{x,t) 



S{x,t) 



(vY 



(145) 



Using the estimate (|49l) proved in lemma [U with m > 3, we obtain a bound for the first 
term in ( 11421) 



4vrf 



0(2/) d(T (y) 



(t+\x\)(t-\x\) m - 2 



>S(x,t) 

for almost all t G R + and all x G R 3 . 

We proceed analogously with the second term in (11421) . with m > 3, using the estimate 
(150]) from lemma [1] for the term containing /, and obtain 



1 

4^2 



(y - x) ■ Vf(y) da(y) 



S(x,t) 



< 



1 

4"vrt 



\Vf(y)\da(v)<fi 



S(x,t) 



r 0) 

(t+ \x\)(t- \x\) m ~ 2 

(147) 



6 By "variation of Fubini's theorem" we mean here the non-standard decomposition of IR 3 = S(x, ■) x R + . It 
is equivalent to introduction oig x (u,t) := g(x + ujt) defined on 5(0, 1) x R + , where the standard Fubini's 
theorem holds. 



22 



1 



4vrt 2 



f(y) da(y) 

'S{x,t) 

for almost all t G and all iGl 3 . The last three estimates lead to 



~ f °(t+\x\)i-\x\)m-2 ^ 



\v(t x)\ < Cm -^o + fi + fo) (U9) 

with C m defined in (11401) . for almost all t G M+ and all x G R 3 . The inclusion of the 
boundary t = is trivial, because v(0, x) = 0, although not necessary, since the set {0} x M 3 

1+3 

is of measure zero in M + . Thus 

IK^IU^ = \\(t+\x\)(t-\x\) m - 2 v(t,x)\\ Loc(K+3) < C m . (150) 

It remains to show that v is measurable on R + . A product of two measurable functions 
is measurable, hence g(y) as well as (y — x)V/ and / are measurable on M 3 . The solution 
v(t, x), given by (I142p . can be rewritten as 

f tg{x + tu) + tu ■ Vfjx + tu) + f{x + tu) 
v(t,x)= da{u) (151) 

JS(0,1) 47r 

where we have changed variables so that we integrate over a unit sphere around the origin. 
Define h(t, x,u) := [tg(x + tu) + tu ■ V/(x + tu) + f(x + tLj)]/Air, which is measurable on 
K + x IR 3 x S(0, 1), because it involves only sums, products and compositions of measurable 
functions. Then, 



v(t, x) = / h(t,x,u) da(u), (152) 
Js(o,i) 

according to the theorem of Tonelli, is a measurable function of (t, x) G M + x M 3 defined 
almost everywhere. Together with boundedness almost everywhere, shown in (11501) . we 
conclude that v G Lf m _ x . 

Since v is defined by (11421) almost everywhere and v G L lrjc , by lemma [31 v solves weakly 
the free wave equation (|6T|) . □ 



Lemma I5al 

Proof. v(t,x), given for t > 0, x G M 3 by the formula 

v{t,x):=^ [ ^^-d K {s,y) (153) 

]__l_3 1+3 

for F G C°(IR + ), is well-defined everywhere in (t, x) G M + x IR 3 = R + , because the integral 

1+3 

over a compact set if t) G R + of a continuous integrand exists and is finite. Then, using 
a variation of Fubini's theorem, we split the integration 

v(t,x) := — / ds ( da(y) 



4Wo JW„ t- S (i54) 



— [ dt f da{uj) (t-s)F(s,x + (t-s)uj). 

4?r io 75(0,1) 
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Now, H(t, x, s, uj) := {t — s)F(s, x + (t — s)u>) / (47r) is continuous in all variables t, x, s, ui and 
for all (s, uj) G [0, T] x S{0, 1) is uniformly bounded on (t, x) G [0, T] x E 3 by \H(t, x, s,u)\< 
TF G L 1 ([0,T] x S^O, 1)) for any T > 0. Then, by a standard theorem on Lebesgue 
integration of continuous functions, the integral 

v(t,x)= ds da(ui) H(t,x, s,lu) (155) 

Jo JS(0,1) 

1+3 

is continuous on (t, x) G [0,T] x M 3 for any T > and hence on (t,x) G 1R + . (The case 
t = is trivially included, because lim t _ >0 + v (ty x) — — v(0, x).) 

Since v is defined everywhere in M+ +3 and hence almost everywhere in R + by the above 
formula, and v G L\ oc , then by lemma [3] it solves weakly the wave equation (1641 ). 

Now, knowing that t>(t,x) is continuous and hence finite everywhere, we can estimate 

Ht,x)\ < J- / S^1^( S)2/ ) 

47T /i<7 r /\ I — S 

< F [ dK{s,y) 



4?r y^^t) (j/>«(s + - \y\Y~\t - s) 

Lemma [2] gives an estimate for this double integral (integral over a cone) and we obtain a 
bound 

\v(t,x)\<- , C ff F ° , x : (157) 

which is valid for all (t,x) G M. 1 ^ 3 . This bound is also satisfied for t — 0, because f(0,x) = 
0. □ 



Lemma _ 

14-3 

Proof. We have G L^ p then also F G L°°(M + ) and hence F G Li oc (R + ). The 

14-3 

whole space-time R + can be covered by the set of cones {K(x,t) : t G up to one 

point (0,x) G M_)_ x M 3 , which itself plays no role in integration and measurability, i.e. 

jig 

W + = C(x, •) x M + . Then, from a variation of Fubini's theorem 7 it follows that F\x(x,t) G 

L 1 (-ft'(x, t)) (we skip "local" because K(x, t) is compact in M^ 3 ) for almost all t G R+ and 
all x G M 3 . Then, of course, F(s,y) and hence F(s,y)/(t — s) is measurable on the same 
cones, so the following integral over the cone 

v(t,x):=^-f ^AdK{s,y) (158) 



4?r .lK{x,t) t — s 



exists for almost all t G IR4. and all x G 



7 By "variation of Fubini's theorem" we mean here the non-standard decomposition of R 1 ^ 3 = C(x, •) x M + . 
It is equivalent to introduction of F t . x (uj, s) := F(t—s, x+uis)) defined on S(0, 1) x [0, t], where the standard 
Fubini's theorem holds. 
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As next, its value can be estimated. The condition \\{x) q F(t,x)\\Lf p = Fq < oo means 
that \F(t,x)\ < F /((x) q (t + \x\)(t - |z|) p_1 ) for all (t,x) G !+ +3 except a set A of measure 
zero in IR + . From a variation of Fubini's theorem 7 for measurable sets and the decompo- 

— 1+3 

sition ]R + = K(x, •) x R + it follows that for almost all t G M. + the measure (on the cone 
K(x,t)) of A\k{ x ,£) is zero. It means that for almost all t G M+ and all x G M 3 it holds 
\F{s,y)\ < F / ((y) q (s + \y\)(s— almost everywhere on the cone (s,y) G K(x,t). For 

these (t, x) we can estimate 

Mt,x)\<±- [ l -^fdK(s,y) 

47T jK{Tf\ \t — S\ 

(159) 



< 



'(*,*) 

F /" dK(s,y) 



4vr 7 A>jt) (j/)9(a + |j/|)<a - - s) 

Lemma [2] gives an estimate for this integral over the cone and we obtain a bound 

\v(t,x)\<- / C ™' F ° „ : (160) 

(t+ \x\)(t - \x\yp- 1 v ; 

which is valid for almost all t G M+ and all x G M 3 . (This bound is also trivially valid for 
t — 0, because w (0, x) = 0.) Finally, it gives 

< C M F . (161) 

In order to have t> G L^° p it remains to show that v(t, x) is measurable on IR + . To this end 
we use again a variation of Fubini's theorem with K(x,t) = S(x,t — •) x M + for (11581) and 
those (t, x) where this integral exists and is finite, and write 

v (t,x) = —[ dt f da(u) (t- s)F(s,x + (t- s)w) (162) 
47r Jo Js(p,i) 

Since (x) q F(t, x) G the function F(t, x) is measurable on M + x K 3 . Then H(t, x, s, u) := 
(t — s)F(s, x + (t — s)uj)/(Ait) is measurable on (t, x, s, w) 6 R + x I 3 x M + x ^(0, 1). From 
Tonelli's theorem we get 

v(t, x) — dt da(u) H(t, x, s,lu) (163) 

Jo Js(o,i) 

1+3 

measurable on M + x R 3 = R , . Finally, v G L\ oc and is defined almost everywhere by (11581) , 
hence by Lemma [3] v solves weakly the wave equation (I67|) . □ 



Theorem | 

Proof. Consider the following iteration scheme 

u n+1 :=I (f,g) - L (Vu n ), n = 0,1,2,... and u_i := 0. (164) 
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For g,Vf G and / G w ^ n m > 3 from lemma lib! we get u = I (f,g) G ^i° m _i. 

Next, observe that if w n G with some p > 1 then 

||(x) fc yU n || L oo ) < IK^Vlli^ll^Hio^ = VoHUnlU^ < OO (165) 



and from lemma [5b] with F = Vu n we get Lo(Vu n ) G Lf when p < k. Because Lf pi C L° 
when pi > p 2 , we get u n+1 G with p < min(m — 1, k). By induction we obtain u n G Lf p 
for every n = 0, 1, 2, ... with the optimal value p := min(m — 1, k). Then, we have 



\u n+1 - u n \\ L ^ p = \\L (-V(u n - u n _i))||z°° p < C P)k \\{x) k V[u n - u n _i)|| L - 
< C P: k\\ (x) k V\\L°c \\u n — u n -x\\i,oo p = C P: kV \\u n — u n ^i\\ L ^ 



again making use of lemma l5bl with F = —V{u n —u n -\) G (x)~ h L^ p . For 5 := Cp^Vo < 1 the 
iteration is a contraction in Banach space Lf p . A simple argument shows that the sequence 
u n is Cauchy. We have 



IK+i - u fc || L » < <5 fc+1 ||u - u-xWl^ = 5 k+l \\hU\g)\\ Lrp (167) 

and for n > m 

n—m—l n—m—l 
\\u n ~ U m \\ L oo p < ^ \\ u m+k+l ~ U m+k \\ L oo p < Y 5fc+m+1 |l J o(/>^)IU^ 

fc=o fc=o (168) 

Xm+1 

and this expression can be made arbitrarily small (< e) for all n,m > M(e) when 5 < 1. 
Since Lf^ is a Banach space the Cauchy sequence u n has a limit u G with the property 

u = I (f,g)-L (Vu) (169) 

and hence 

< ||/oC/,<?)||l~, + ||^(Vu)|U~ < C m (fo + fi+go) + C p , k V \\u\\ L ~. (170) 

Then 

I, || . Cm(/o + /l + go) _ n n „ 1 \ 

Moreover, equation (11691) together with lemmas Ha] and [5a] imply 

(it / d 3 x Dip(t, x) u(t, x) = - / d 3 x d t <p(0,x) f(x) + / d 3 x <p(0, x) g(x) 



(172) 

(it / d 3 x (p(t, x) V(x) u{t,x), 



what gives equation (Q. Since u G L[ oc it is a weak solution to the wave equation ([I]). 
Uniqueness is guaranteed by theorem [[] □ 



26 



Theorem I2al 



Proof. The condition |u(t,x)| < C/((t + \x\)(t + |a;|) p l ) implies u G Lf p and together 

with u G C°(M + ) these two conditions are equivalent. Therefore, proof of this theorem is 
analogous to the proof of theorem 12 bl with the only difference that continuity must be shown 
separately in the iteration. 

Analogously, for iteration defined by 

u n+ i := Io{f,g) - L (Vu n ), n = 0,1,2,... and u_i := 0. (173) 

and (f,g) G C 1 ^ 3 ) x C°(M 3 ) satisfying the bounds (J3TJ, hence g,Vf G Z£, by lemma HH 

i I o 1-1-3 

we get M = W,g) G L^ m _ 1 nC°(M + ). Next, if u n G L^ p nC°(M + ) for some p > 2, hence 
also Vw™ G C°(IL 1 ( +3 ), because V G C (M 3 ), then analogously, by lemma EH with F = Vu n , 
we get L (Vu n ) G Lf p n C°(I+ +3 ) when p < fc. It implies that u n+l G L~ n C°(I+ +3 ) with 

p < min(m — 1, k). By induction we obtain u n G L^° p nC°(]R^ +3 ) for every n — 0, 1, 2, ... with 
the optimal value p := min(m — 1, fc). Analogously, using lemma [Sal with F = V(u n —u n _i) G 
(x)~ k L™ p n C°(I+ +3 ) we show 

||«n+l - M n||L- p < C pik V \\u n - U n -l||i^ (174) 

1-1-3 

what presents a contraction in Banach space Lf 3 fl C°(K + ). Again by full analogy, for Vq 

1-1-3 

small enough, we show that the sequence u n is Cauchy in Banach space Lf 3 fl C°(R + ), 

1+3 

hence u n has a limit u G Lj° p fl C° (R + ) with the property 

u = I Q (f,g)-L (Vu). (175) 

Then 

1+3 

Since u G C°(1R + ) the last condition implies 

\u(t,x)\<-. , , W C , „ : V(t,x)GM+ +3 . (177) 

And also analogously, equation (11751) together with lemmas [4b] and [5b] imply 

(it / cftr n<p(t, x) u(t, x) = - / d 3 x d t (p(0,x) f(x) + / c/ 3 x y?(0, x) g(x) 



(178) 

c/t / d 3 £ (p(t, x) V{x) u(t,x), 



what gives equation (Q. Since m G L[ oc it is a weak solution to the wave equation ([I]). 
Uniqueness is guaranteed by theorem [0 □ 
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Corollary Q] 



Proof. Analogously as in proofs of the theorems above we define an iteration 

u n+ i := I (f,g) -L (Vu n ) - L (F), n = 0,1,2,... and u- X := 0. (179) 

From lemma HE] we get u = I (f,g) £ Lf > m _ 1 . From lemma [5b] we get L (F) £ L^ r . 
By same argument as in previous theorems we have Lo(Vu n ) £ Lf> with p < k provided 
u n £ with some p > 1. Then, by induction we obtain u n £ L^ p for every n = 0, 1, 2, ... 
with the optimal value p := min(m — 1, k, r). Then, we find 

\\u n+ i - u n \\ L oo p = \\L (-V(u n - tt n _i))||z~ p < C , Pifc ||(x) fc V(M n - u n -i)\\ L ^ 

k (^"^) 
< Cp,fc||(^) ^||l°°||w„ — Wn-l|Uf° = Cp,fcVo||Un — u n-l||Lf° p 

again making use of lemma l5bl with F = —V{u n —u n -i) £ (x)~ k Lf p . For 5 := Cp^Vo < 1 the 
iteration is a contraction in Banach space and the sequence u n is Cauchy. It converges 
to the limit u £ (Banach space) which satisfies 

u = I (f,g)-L Q (Vu)-L (F) (181) 



and hence 



|«IU- < ||/ (/ ) ^)||l- + HMVVlll^ + ||L (F)|| L oc 
< C p (/o + fx + 9o) + C P!k V \\u\\ LT + C m F Q . 



'1821 



Then 



IHI ^ ^ = a (183) 

To prove the continuous case one needs only to use the lemmas Ha] and [5a] instead of the 
lemmas [4b] and [5b] and repeat the convergence argument from the proof of theorem [2al □ 



VI. DECAY ESTIMATES FOR THE DERIVATIVES 

One can state the estimate given in theorem [2a] or [2b] in a more detailed form, namely 
that the solution of ([I]) satisfies 

C m - , + ||V/||loc + II^IUoo) 
||ti|| L oc < v " " " " m " " m> (184) 



t-c p , k -\\v u 



k 



when k > 2, m > 3, p := mm(k,m — 1), the norms H/Hz, 00 ^, ||V/||l^, ||<7||.l°j are finite and 
||V|U- < C-J < oo. 

If u is a classical solution one can differentiate the wave equation (TTJ with respect to time 
and obtain 

Ud t u + Vd t u = 0. (185) 

If u is a weak solution of ([ID then this equation is also satisfied in the weak sense. The 
existence of the weak solution will be guaranteed by the theorems below. It can be treated 
again as equation ([rj for a new variable v := d t u with the initial data 

w (0, x) = g(x), d t v(0, x) = d$u{0, x) = -Af, (186) 
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where A := —A + V(x). Then, theorem [2a] or [2b] gives existence of the solution dtu G 



C°(IR + ) or d t u G respectively, and the estimate 



PHIl- < : ~ nTTTi ■ ( 187 ) 



1-C^-||^|| L 



fe 



The term 11^4/Hz, 00 , can be bounded by ||A/|| L ^ + ||V/||.l°°,- Analogously, one can obtain 
existence and an estimate for d^u 

C m -{\\Af\\ LZ + ||V(Af)|U~ + \\A g \\ L -,) t 
Uu\\l Tp „ < fluF — ■ ( 188 ) 



The term ||V(^/)|U- can be bounded by ||VA/|| i% + ||^V/|U~ + ||(W)/IU°v The 
three estimates can be put together when the constants m,m',m" and p,p',p" are related 
to each other. Using the following property of the weighted L°° norms \\hi ■ /^IIl 00 ^ < 
||/ii||l«> • H^lUg we get for m = m! = m" > 3 

Corollary 2a. 

u, d t u, d 2 t u G L™ p (189) 

with p := min(/c, m — 1) provided 

/, V/, A/, VA/, V</, G L~ (190) 

||V|| L ~ < C~l < oo, llVyllx- < oo. (191) 

Alternatively, choosing m = mf — 1 = m" — 2 > 3 we get for classical solutions especially 
simple form of the assumptions 

Corollary 2a'. 

IMU&j \\dtu\\ L ™ p/ , \\dtu\\ L oc p// < oo, (192) 

with p := min(/c, m — 1), p' := min(/c,m), p" := min(/c, m + 1) provided (f,g) G C 3 (M 3 ) x 
C 2 (M 2 ) m& 

|V 3 /I G ^ +2 , |V 2 s| G L- +2 (193) 

and V G C X (M 3 ) 

||V|U- < min^,^,^) < oo, \\VV\\ LT < oo. (194) 
Here, we introduced a simplified notation 8 

3 

|V7>|:= \ d ai- d «nh\. (195) 

ai,...,a n =l 

Because of the regularity, the bounds on |V 3 /| and |V 2 g| can be integrated to give 

||V/|U~, ||A/|| L oo +i , ||VA/||^ +2 < oo, (196) 



This notation simplifies the assumptions, but looses information about the decay of directional derivatives, 
which is however rarely used. 
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and 

Nk«, l|v^|U- +1 , \\&g\\L% +2 < oo. (197) 

The control of spatial derivatives is more difficult, since V depends on x. Again, by 
differentiation of the wave equation ([T|) with respect to x l we get (for weak solutions the 
equation holds as well by the same argument as above) 

UdiU + VdiU = -(diV)u. (198) 

This equation can be handled with Corollary [1] for v = diU and F = —(diV)u with the initial 
data v(0,x) = dif(x) and d t v(0,x) = dig{x). Then, 

CW(||<9i/|U~, +||V^/|| L ^ + ||^||^ / ) + C p ^||(x) fc '(^y)n|| Lr 
II^mIIl-, < : m „ llT , M ~ (199) 



k 



for p' := min(/c, vnl — l,p) and m' > 3, k, k' > 2, k' > p > 1. The last term in the nominator 
can be bounded by ||(9iV)||x,^||w||Lj° . Analogously, an estimate for <9jc\u reads 

\\d i d k u\\ L ^ n < ^"(ll^fc/IU",/.! + l|V<9i<9 fe /|| L ^ ;/ + ||9i9fc^||£«J 

for p" := mm(k,m" - l,p,p') and m" > 3,k,k',k" > 2, k" > p > l,k' > p' > 1. We can 
combine these estimates for m — m' — m" > 3 and k = k' = k" > 2 and get 

Corollary 2b. 

it, |Vu|, |V 2 u| G L^ p (201) 

iwtt p := min(/c, m — 1) provided 

f G L%_ x , |V/|, |V 2 /|, |V 3 /U, |V^|, |V 2 ^| G L~ (202) 
V, |W|, |V 2 V| G and \\V\\ Lf ° < C~£ < oo. (203) 

Finally, combining corollary [2a] with [2b] we get corollary [2] 

The same reasoning can be repeated for (/, g) G C 3 (M 3 ) n C 2 (1R 3 ) and 1/ G C 2 (R 3 ) which 

1-1-3 

give the classical solution u G C 2 (R + ). The only difference is that one needs to use theorem 
a] instead of [2b] and the continuous version of corollary [TJ 
Higher derivatives can be treated in an analogous way. 
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